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Abstract 

A Wigner function representation of multi-band quantum transport theory is devel- 
oped in this paper. The equations are derived using non-equilibrium Green's func- 
tion formulation with the generalized Kadanoff-Baym ansatz and the multi-band 
k.p Hamiltonian including spin. The results are applied to a two-band resonant 
inter-band tunneling structure. 
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1 Introduction 



The single band approximation [2] is the most often used approach in quantum 
device models. In this approximation, inter-band processes in the structure are 
ignored and the boundary conditions for the model are oversimplified. The 
single band electron transport models have been applied to large band-gap 
semiconductor heterostructures (e.g. AlAs/GaAs/AlAs). 

Multi-band quantum transport has attracted attention due to the existence 
of various inter-band tunneling structures and especially resonant inter-band 
tunneling structures (RITS). It is possible to achieve multiple negative differ- 
ential regions (NDR) in these structures and they are shown to exhibit high 
peak-to- valley current voltage characteristics. RITS's are based upon the type- 
I, the type-II staggered and the type-II broken-gap band alignments. Recently 
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[3], it was theoretically shown that a type-II staggered band-gap resonant tun- 
neling diode can exhibit oscillations in the THz region. A significant amount 
of inter-band current can be present in a staggered band-gap structure. The 
coupling between the conduction and the valence bands is considered to be 
the dominant mechanism in these structures. Therefore, in a correct descrip- 
tion of electron transport, multi-band effects must be included [4]. We do not 
consider an external magnetic field but the extension is quite straightforward. 

For direct band-gap semiconductors the conduction band near k = has 
symmetry properties (spherical symmetry) same as the \S > atomic orbital 
(/ = 0,mi = 0). On the other hand, the valence band near k = has sym- 
metry of p-orbitals, \X >, \Y >and \Z > (p-orbitals are antisymmetric and 
I — l,mi — —1,0,1). An eight-band model [5] can be gained with the in- 
clusion of spin. These states that become doubly degenerate, \S T>, \X |> 
, \Y t>, \Z t>, \S i>, \X J>, \Y i>, \Z i>. The spin-orbit interaction lifts the 
six-fold degeneracy of the valence band and splits it into a four-fold degen- 
erate and a two-fold degenerate level. If the spin-orbit coupling is considered 
in the energy band calculation of a semiconductor, the Bloch states become a 
mixture of spin up and spin down states. This becomes important in asymmet- 
ric quantum well devices. Inversion asymmetry of the bulk or the confining 
potential causes spin splitting even in zero magnetic field due to spin-orbit 
interaction. Therefore inclusion of spin and spin-orbit interaction in the quan- 
tum transport equations is important. 

Derivations of quantum transport equations have usually been based on the 
first-order gradient expansion [6] . This approximation is based on the assump- 
tion that the "fast" quantum variations can be separated from the "slow" 
macroscopic variations and causes the loss of information related to quantum 
processes such as interference and tunneling which are crucial in nano-scale 
devices. Buot and Jensen [2], [7] presented an alternative derivation for single- 
band devices and provided an exact integral form of the quantum transport 
equation which is capable of accurately describing full quantum effects. The 
first-order gradient expansion is still needed to simplify the collision terms af- 
ter the derivation of transport equations. Their approach has been generalized 
for the multi-band transport in this work. 

The Wigner function modeling of the quantum transport in the single-band 
resonant tunneling structures has been quite popular in the literature due to 
its success in dealing with the dissipation and the open boundary conditions 
[2]. Similarly, it is expected that one should be able to model the multi-band 
transport in the resonant tunneling structures using the Wigner function. 

There are a number of works on multi-band Wigner function representation of 
quantum transport in the literature. Miller and Neikirk [8] used Wigner func- 
tion for multi-subband transport in double barrier resonant tunneling struc- 
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tures. Demei et al. presented multi-band Wigner function formulation without 
spin [9], [10]. Zhao et al. [11] showed that multi-band quantum transport equa- 
tions can be decoupled to reduce the number equations to be solved. Borgioli 
[12] employed a two-band Kane model to derive Wigner function equations 
for resonant inter-band tunneling diodes. 

The point of this work is to develop a complete theory of the multi-band 
Wigner function for transport in nano-scale devices. This has been accom- 
plished by using non-equilibrium Green's function methodology which is known 
to be the most complete description of quantum transport. The results give us 
the Wigner function formulation of multi-band systems based on k.p theory. 
The results can be easily simplified by the symmetry arguments of the band 
structure of the system under study. The derived multi-band Wigner function 
equations which are also capable of description of zero magnetic field spin 
transport devices are the first in the literature. 

In the introduction, the subsection 1.1 we present a preliminary on k.p method. 
Then in the subsection 1.2, the non-equilibrium Green's function method in 
phase space is given. We derive the Wigner function equations for multi-band 
systems in section 2. Finally, in part section 3, we apply the formalism to a 
simple one dimensional two-band resonant inter-band tunneling diode. 



1.1 k.p Hamiltonian 



The Schrodinger equation for the lattice periodic part of the Bloch functions 
can be written as [13] 



+ V(j) + 7\ + — k.p+ -^-(a x v^)-P 



2m 2m m Am^c 1 

+ \ 2 (o- x vV).k]|nk >= £„ k |™k > . 



We express the bulk band matrix element H ab of the Hamiltonian in the second 
order k.p theory as 

H ab = D^k a k, + D$*k a + (D^ + Va (v))5 ab (2) 



where the indices a and 7 are summed over x, y, and z. The o, b include 
both the band and the spin indices. V(r) is a spin-independent self-consistent 
potential. Note that for heterostructures k a is replaced by — i\J a . 
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We define a vector n as, 



So, 



(2)q 7 
a6 



2m n 



^ab^a'y "I - 



'A 

v m 



) 2 E 



a 7 



E a +E b 



Er)' 



(4) 



noting that the second term arises from Lowdin renormalization and needed 
to include the interactions with remote bands. So we denote these states by 
the index r. These interactions are usually ignored in Kane models and so 
D^ ai {a 7^ b) terms vanish. 

The part of H a i, linear in k include the inter-band coupling (k.p interaction) 
and the spin-orbit interaction terms, 



where 



Kb =< U a \p a \U b > +^—; < U a \{a x V V) a \U b > (6) 



for a 7^ b. Note that 7T ba = (n^ b )* and 7r" a = (which implies that D^ Q terms 
vanish) . 

The terms p^ b =< U a \p a \Ub > are the inter-band momentum matrix elements 
and measure the strength of the coupling between the various bands. Note 
that 7r" a = even if the band minimum is at some point other than k = 0. 
The term ■^ j V a \ * s usually written in terms of a real parameter P originally 
defined by Kane. The value of this parameter is known for any given material. 

in 1T1 ?f) 
P = < S\p x \X >= < S\p y \Y >= < S\p z \Z > . (7) 



The band edge is denoted by Df® such that, D a °J = E a (k = 0). We write 
D a °a + K(r) as E a (r) in the calculations. 
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1.2 The non- equilibrium Green's function formalism in phase space 

The multi-band Green's function [1] is defined by 

G a6 (l,2) = -i<^ a (l)^ t (2)> (8) 

where C denotes that time arguments are on a contour rather than a real- 
time axis. The expectation value is defined in a grand-canonical ensemble. 
We define the space time arguments 1 = (ri,ti), 2 = (r2,t2). ip a is the field 
operator for electrons. The a, b include both the band and the spin indices. 
The equation of motion of the band electron Green's function is written as 
(sum over repeated indices) 

\ih5 a p—-H aP {l)]Gp b {l,2) = 8 ah 5{l-2) + | d3£ a/3 (l, 3)6^(3, 2), (9) 

and the adjoint equation is given by 

HM a/3 — - H a p(2)]Gp b (l, 2) = 5 ab 5(l - 2) + J d3G af3 (l, 3)^(3, 2), (10) 

where the self-energy is denoted by £(1,2). It describes the scattering of 
electrons by other electrons, phonons and impurities. Throughout the paper, 
Greek indices are used to denote the repeated indices to be summed over. 

The generalized Kadanoff-Baym (GKB) equation describes the time evolu- 
tion of the electron correlation function G< a (l,2) in the band a. It should be 
noted that, in the equal time limit (ti = t 2 ), the off-diagonal (a ^ b) lesser 
Green's functions correspond to inter-band polarizations in energy band space 
and inter-spin-band polarizations in spin space whereas the diagonal Green's 
functions give the particle densities with the spin up or down in each band. 
The time evolution of G < (1,2) given by the GKB equation can be written, 
using the Langreth algebra, as 

m{ ~k + W 2 )G « b{1 ' 2) = [H ' G<](1 ' 2)ab + [s< ' ReG * ](1 ' 2)ab 

+^{£<,A}(l,2) o6 -l{r,G<}(l,2) a6 . (11) 

where [ ] is the commutation and {} is the anti-commutation. Above, the 
spectral function is defined as, 

4*(1,2) = <[G^(1,2)-G<(1,2)] = -2ImK(l,2)], (12) 
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and the dissipation function is 

r ab (l,2) = i[E> - £<] = -2Im[£*(l,2)]. (13) 

Switching to center of mass and relative coordinates (Wigner coordinates) 
done by defining 



ri + r 2 _ h. +t 2 

2 ' ~ 2 

v = r 2 - ri; i = t 2 - h- (15) 



R=£l±£i ; r=^, (14) 



The four-dimensional, (3+1), crystal momentum and its conjugate variable 
lattice coordinate are represented as p = (p,E), r = (R,T). Note that we 
use hh as the crystal momentum when the matrix elements of the Hamilto- 
nian is considered. Therefore Ti\t and p are used as the crystal momentum 
interchangeably. 

The Weyl- Wigner representation, W[0] = 0(p,r), of any operator 0(1,2) 
[14], [15] is defined by 

0(p,r)= I dvexp(^p.v) < R - ^-\0\R+ ^ > . (16) 



It is very important to note that 0(p, r) is real if O is hermitian. Let C = AB 
then the differential form of the Weyl transform of the product of two operators 
can be obtained by, 

W[C] = C(p, r) = exp(ik)A(p, r)B(p, r) 

= eM-iA)B(p,r)A(p,r), (17) 



where 



2 dr dp dp Or 



The partial differential acts on only A and acts on B only. 

To obtain the Wigner function equation, it is necessary to switch to a phase 
space description of GKB. Taking the Weyl- Wigner transform of both sides 
of the equation (11) gives the GKB in the phase-space-energy-time domain 



i^G<(p, r) = exp(*A) [H, G<] (p, r) ab + exp(*A) [£< , ReG R ] (p, r) ab 
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+ l - exp(*A){£<, A}(p, r) ah - % - exp(*A){r, G<}(p, r)09) 

For any operator A and B, the integral representations of exp(iA)A(p, r)B(p, r) 
and exp(— iA)A(p, r)B(p, r) in (3+1) dimensions can be written as [1], 

i f i i 

exp(±iA)A(p, r)B(p, r) = jj^y J dridpidr 2 dp 2 exp[-p l .(r - r 2 )] exp[-ri.(p - p 2 )\ 

xA(p±^,r T ^)B(p 2 ,r 2 ), (20) 

Defining 

dp x dr x exp(-pi.(r - r 2 )) exp(-r 1 .(p - p 2 ))A(p ± — , r =f — ), (21) 

the equation (19) becomes 

d 1 f 

ih 0fG< b (p, r) = J dp 2 dr 2 K c H ^(p, r-r 2 ;r,p- p 2 )G^ b (p 2 , r 2 ) 

+ ~{hf I d P2 dr 2 K x< f} (P> r ~ r 2',r,p - p 2 )ReG^ b (p 2 ,r 2 ) 
+ 2{j^y I d ^ driiK yj<P^ ~ r 2;q,p- P2)A pb (p 2 ,r 2 ) 
~2(h^f I d P2dr 2 K s T ^{p,r - r 2 -r,p - p 2 )Gj b {p 2 ,r 2 )i22) 

where K s /(p,r-r 2 ;r,p-p 2 ) = K^(p,r-r 2 ;r,p-p 2 )±K^(p,r-r 2 ;r,p-p 2 ). 

The multi-band Wigner function is found by taking the energy integral of the 
Weyl-Wigner transformed G^: 

/ a6 (p,R,T) = J dE(-i)G<(p,E;R,T). (23) 

Note that the indices a, b include both spin and band. The total Wigner 
function of the multi-band system with spin can be written as the summation 
over the band and the spin indices, 

/(p,r,«) = X; £ <^/ c r'(P,r), (24) 

c,d m,m'=1,l 

where c, d are band, m and m are spin indices, k takes values of x, y, and z. 
a is the unit matrix and the others are the Pauli spin matrices [16]. Therefore 
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each intra-band and inter-band component of Wigner function becomes 2x2 
matrix in spin space. 



2 The Wigner Function Equations for Multi-band Systems 



Under the assumption that the self-energies are slowly varying with respect 
to the center of mass coordinates, equation (22) reduces to 



d 1 r 

in Qf G ab(P> r ) = Jj^y J dp 2 dr 2 K c Haf3 (p, r-r 2 ;r,p- p 2 )G^ b (p 2 , r 2 ) 

+E> /3 (p, r)G% (p, r) - £< (p, r)G> b (p, r) (25) 



using 



The self-energy function can be written as [17], 

2^(1,2) =iG>'<(l > 2)D>'<(l J 2). (27) 

The Weyl- Wigner transform gives of the above equation (27) gives, 

E^for) = dqG>f{p + q)D>><{q). (28) 

Assuming the phonon bath is in equilibrium, the Fourier transforms of the 
phonon Green's functions can be written as, 

D<(q, E') = -i/iMq[(iV q + l)8(E' - Q q ) + N^E' + ft q )], (29) 

D > ( q , e') = -ihAQ(N^ + i)S{e' + n q ) + iv q <5(£' - n q )] (30) 



where M q is the electron-phonon scattering matrix element. Therefore, inclu- 
sion of the phonon scattering gives the following scattering functions 

^ab= E i ^ / dqG^(p + q, E + ^ q , r)JUj(J\T q + 1 + ± v ), (31) 
S> 6 = _E i 7^ / ^> b (p + q, £ + r)Mq(Nq + (32) 
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The first term on the right hand side of the equation (25) can be written as 

exp(iA)[H, G < ](p, r) ab = exp(iA)H a/3 (p, r)G| b (p, r) - exp(-iA)Hp b (p, r)G^(p, r).(33) 

The integral representation, using the equation (20) becomes 

If % % 

exp(iA)[H,G < }(p,r) ab = J dr 1 dp 1 dr 2 dp 2 exp( y -p 1 .(r - r 2 )) exp(-ri.(p - p 2 )) 

x[H aP (p+^, q -^)G< b (p 2 ,q 2 ) 



-Hp b (p-^,T+ r -±)G<p(p 2 ,r 2 )}, (34) 



where 



+D ( »°(p a ± P f) + (D<»> + K(r T ^))^. (35) 



At this point, since the purpose of the paper is to derive a Boltzmann type 
transport equation, it is useful to make quasiparticle approximation to get the 
form of the spectral function. The free generalized Kadanoff-Baym (FGKB) 
ansatz for multi-band systems is stated as [18], 

- iG< (p, E, r, T) = hf ab (p, r, T)5(E - ^M±EM) (36 ) 
iG> b (p, E, r, T) = h(5 ab - Up, r, T))S(E - EM±^M ) (37) 



Using FGKB ansatz, one can simplify the scattering functions so that the 
equation of motion for G < in the phase-space-energy-time domain becomes 



ih-^G<(p, E, r, T) = D^ PaPl G< - D^p aPl G< 



gFp^abyVi 1 > ± J — ^af3 VaFy^ph ^ (] b W7 U a/3 



ft 2 n (2) a7 ^ a < (2)Q7 a a < 
"T^ a^a^ + T 1 ^ a^a^ 

+»a? \Pa + 2idr J<*l* D p b \P* ^arj^ 

+<W^ / dP2dri exp(^n.(p - p 2 ))[Df} + Ka(r - y)]G< (p 2) r) 
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~ 5pb h I dp2driex v^ ri -( p ~ p2 ^ D w + y w( r+ 

+ Pa/3^|b ~ S a/3^] (38) 

where p a denotes p x ,P y ,Pz , -£r denotes J^, J^, ^ and there is summation 
over Greek indices as usual. Using the equation (23), the multi-band Wigner 
functions for an open system in weakly contact with a phonon heat bath can 
be written as 



d 

ih—f ah (p, r, T) = D^pap^fatp, r, T) - Df^p^f^p, r, T) 



h 2 _f 2 w d d h 2 ( 2 ) Q7 a a 

-^(p,r,T) + T DV — - 



+ |^]^(P>r,T) - D$ a \p a - l£-] faf) (j>,T,T) 
y rfp 2 ^i exp(-ri.(p - p 2 ))[ J Di° ) + Ka(r - y)]//3b(p, r, T) 

~ 6pb h I dp2driexv ^ ri -( p ~ p2 ^ D w 

+ Z ^/ rf q(^/3 - /a/3(p + q,r,T))/ /35 (p,r,T) 
J ( l'(p + q) + ^(p + Q) _ S,(p) + *( P ) + + 1 _ 1 

- £ ^/rfq/«/3(P + q.r,r)(«»-/^(p,r,T)) (3 
J ( l(p + q) + ^(p + q) _ S,(p) + *( P ) + ^ q)M , (jVq + 1 + 1 

The third and the fourth terms at the right of the equation (39) are the usual 
drift terms. The first, the second, the fifth and the sixth terms do not cancel 
each other for o ^ b if and only if Lowdin renormalization is considered. If the 
effects of the remote bands are ignored, these terms cancel each other. The 
seventh and eight terms explicitly give the k.p and spin-orbit interactions. The 
ninth and tenth terms give the potential term. The last two terms correspond 
to electron-phonon scattering. The relaxation time approximation can be made 
for these. An important simplification occurs when the in-plane wave vector is 
taken to be zero. This approximation gives a set of spin-independent quantum 
transport equations. 
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If the structure under consideration has inversion symmetry as in diamond 
structures, the equations can be simplified further. Note that there is no in- 
version symmetry in zinc-blende structures (bulk inversion asymmetry) and 
spin degeneracy in zinc-blende type heterostructures is lifted even at zero mag- 
netic field. Usually, this splitting is very small and can be ignored. However, 
recently resonant intra-band and inter-band spin filter was proposed based 
on the zero magnetic field spin splitting of the conduction band in the case 
of structural inversion asymmetry (Rashba effect) [19], [20]. We are going to 
discuss the Wigner function modeling of these kind of devices in future papers. 

The total current density is the sum of intra-band and inter-band components 
[14], 

Jir) = iJ dp ^f f ^ r) - (40) 



The particle density in each band is written as 



Let's consider a two-band model (a, b = 1,2) without scattering and neglect 
the effects of the remote bands. Then the first component of the Wigner equa- 
tion (39) becomes 



lh ptt = ~ D ii >«7S-+P 7 7i-]/ii(p,r) 



dT i dr 7 dr a 



If i 
— J dp 2 dr 1 exp(-r 1 .(p - P2)) 



#.('-y)-l'.(' + y)|/„(p !l r). (42) 



The rest of the equations are as follows 



ih df j P ' r) 

or 



2i dr n 



2i dr r 



+ 



If i 
— J dp 2 dr 1 exp(-r 1 .(p - P2)) 



x[^(r-^)-y 2 (r + ^)]/ 12 (p 2 ,r), 



(43) 
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i? .df 22 (p,r,T) _h 



d 



dT 



d 

22 +P7^r]/2 2 (p,r) 



2i dr, 



-]/i 2 (p,r)-L>}r>a- 



2i <9r r 



■]/2i(p,r) 



+ 



/i 3 



y dp 2 driexp(^r 1 .(p - p 2 )) 



x[^(r-|)-y 2 (r + ^)]/ 22 (p 2 ,r) ; 



(44) 



and / 21 = . 



3 1-Dimensional Two-band Kane Model of Resonant Inter-band 
Tunneling Diode 



Resonant inter-band tunneling structures (RITS) are based on the interac- 
tion between the conduction and valence bands and the transport is in the 
growth direction. For narrow band-gap RIT structures, the coupling of the 
in-plane momentum to the transverse momentum component becomes im- 
portant. Therefore a realistic modeling of these structures requires a serious 
amount of computational work. 

The simplest choice is a two-band model that is suitable for large and mid- 
band-gap Type I RITS [21]. The k vector is taken in the z direction and 
the inversion asymmetry is neglected. Therefore, the in-plane momentum 
(k x — k y — 0) vanishes so that the heavy-hole state is decoupled and the 
Hamiltonian matrix becomes block-diagonal [22]. The remote bands are ne- 
glected too. Therefore the Hamiltonian is reduced to a spin-independent three- 
band model (conduction, light and split-off bands) [23] : 



%^p z E lh (z) + £- Q 



3 mo ' 
3 mo " z 



_ flEcv 
V 3 m 





Pz 







E so (z) 



2mo 



(45) 



where P cv = Instead of ignoring the split-off band, Sirtori et. al. 

[23] presented an improved two-band model (conduction and "effective va- 
lence band") that can be gained through a unitary transformation. The 2x2 
Hamiltonian is 



E c {z) + * 



2mo 



mo 



Pz 



—Pz 

mo 1 ^ 



(46) 
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where E v = 2E ^+ E so j g e ff ec tive valence band edge. Therefore the components 
of the Wigner function become 

^ dfcc{Pz,z,t) = %hp z dfcc(p z ,z,t) 
dt m dz 

+ ~hJ rf2; i^2exp(^i(p z -p z2 ))[E c (z- y) - E c (z + ^)}f cc (p z2 , z,t) 
+ZlP f ( n 7f \_J?Lp d fvc(Pz,z,t) 

~T 1 cvjvcyj/zi 1 CV <-> 

m 2m dz 

+^ L p a ,f a ,(p„z,t) + JLp Jf"^*) (47) 

m 2m oz 



.^dfcv(Pz,Z,t) 

lh — at — 



ihpz df cv (p z ,z,t) 



m 



dz 



+ h l dz ± d P^ eX V(-fr Z i(Pz -Pz2))[E c {z - y) - E V (Z + ^-)}fcv(Pz2,Z,t) 



Pz 



ih 



PcvfvviPzi t) i ( .,. 

m 2m 



2 

df m (p z ,z,t) 



dz 



+ P^p f ( n zf \- _!^P 9 fcc(P Z ,Z,t) 

~T -1 CVjCC\fZl n 1 CV r\ 

mo 2mo oz 



(48) 



.^ df vv {p z ,z,t) _ ihp z df vv {p z ,z,t) 
dt m dz 



+ \l dZldPz2eX V(~fr Zl ( Pz -Pz2))[E v (z- y) - E v (z+ ^-)}f vv (p z 2,Z,t) 

-P^P f U Z A + Ap dU(Pz,z,t) 

1 cv J cv \fz j j ~r J cc 

mo 2mo as; 

Z?7 (9 fvci^Pz t Z t^j 

Pcvfvc\Pzi z ,t) — P cv — , (49) 

m 2m Q dz 

and f vc = f* v . Note that above equations are for the conduction and valence 
band electrons and E a (z) = E a (0) + V a (z). The current density can be calcu- 
lated using the equation (40) and given by 



J = Jintra + Jinter = T [ dp z — (f cc (p z ) + fvv{Pz)) + I dpAj ^ — lm[f cv (p z )] .(50) 
n J 777-0 /I J V 2 TYl§ 

where Im[/ C1) ] means the imaginary part of f cv . The particle densities in each 
band using the equation (41) are 

Uc= \l dpz f cc ( Pz ^ ( 51 ) 
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n v = - dp z f vv (p z ). (52) 



4 Conclusions 

In this paper we developed multi-band Wigner function formalism including 
spin, which has a profound effect especially in narrow band-gap semiconduc- 
tors. We employed the k.p Hamiltonian to derive the quantum transport equa- 
tions for multi-band semiconductors using non-equilibrium Green's function 
methodology for systems weakly coupled to a phonon bath. It was shown that 
in addition to drift, potential and scattering terms that exist in single-band 
form of Wigner function, there are terms arising from inter-band coupling 
and spin-orbit interaction. These terms are source of the off-diagonal terms of 
Wigner function in energy band space and spin space. 

A two-band Kane model of resonant inter-band tunneling diode was presented. 
The current and particle densities were derived for this simple model. We are 
going to discuss the numerical solution of the two-band and the three-band 
equations and present the simulation results in future papers. 
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